Current correlations in the interacting Cooper-pair beam-splitter 
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We propose an approach allowing the computation of currents and their correlations in interacting 
multi-terminal mesoscopic systems involving quantum dots coupled to normal and/or superconduct- 
ing leads. The formalism relies on the expression of branching currents and noise crossed correlations 
in terms of one- and two-particle Green's functions for the dots electrons, which are then evaluated 
self-consistently within a conserving approximation. We then apply this to the Cooper-pair beam- 
splitter setup recently proposed^--, which we model as a double quantum dot with weak interactions, 
connected to a superconducting lead and two normal ones. Our method not only enables us to take 
into account a local repulsive interaction on the dots, but also to study its competition with the 
direct tunneling between dots. Our results suggest that even a weak Coulomb repulsion tends to 
favor positive current cross-correlations in the antisymmetric regime (where the dots have opposite 
energies with respect to the superconducting chemical potential). 

PACS numbers: 73.23.-b, 73.63.Kv 74.45.-fc, 72.70.+m 71.10.-w 
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I. INTRODUCTION 

With the development of nanofabrication techniques, 
more experiments in mesoscopic physics draw inspiration 
from quantum optics. Among them, there is a particular 
interest in generating entangled electronic states. Beyond 
the exploration of non-local quantum effects through the 
test of Bell's inequalities in a solid-state device, there 
are potential applications in quantum teleportation and 
information processing, where such sources of entangled 
pairs could be integrated with existing technology and 
infrastructure. 

Because of the spin singlet character of Cooper pairs, 
superconductors constitute a natural source of spin- 
entangled Einstein-Podolsky-Rosen (EPR) pairs. The 
difficulty then lies in spatially separating the constituents 
of these pairs coherently. The commonly proposed setup 
to accomplish such a task is composed of a superconduc- 
tor connected to a fork made out of two normal metal 
leads^"— . The mechanism for such a Cooper-pair beam 
splitter then relies on crossed Andreev reflexion (CAR), 
a process in which the two electrons of a Cooper pair are 
sent into different normal electrodes^"— (as opposed to 
the conventional direct Andreev reflexion which is local). 

Previous attempts on metallic structures revealed the 
difficulty of distinguishing CAR form other contributions 
as there are various parasitic processes^Tiii. Several theo- 
retical works have suggested different directions in order 
to provide preferential enhancement of CAR, from simple 
spin and energy filtering to the effect of strong electron 
interactions, susceptible to favor single-particle over pair 
tunneling ^-SriiS,. 

An important step in the observation of Cooper pair 
splitting came from recent experimental works in tun- 
able double quantum dot devices based on carbon nan- 
otubes^ and InAs nanowires"'^. These systems showed 
great promise as they not only display local Coulomb 
repulsion on the dots, but also allow for the exploration 



of different energy configurations of the dots, two prop- 
erties that should promote CAR over other processes. 

Nevertheless, further efforts are needed both experi- 
mentally and theoretically. Future works should focus on 
time-resolved correlation measurements of the currents 
in the normal leads, particularly relevant to characterize 
non-local effects. However, to this date, they remain to 
be measured experimentally, and still represent an impor- 
tant challenge for theory, as there are only few examples 
of interacting mesoscopic devices driven out-of equilib- 
rium, where current correlations have been obtainedi^Ti^, 
none of which involving hybrid systems. 

In this work, we propose a general approach for com- 
puting currents and current-current correlations in an 
interacting multi-terminal mesoscopic device involving 
quantum dots coupled to multiple normal and supercon- 
ducting leads. While it is knowit2^ that these quantities 
can be expressed in terms of the single and two-particle 
Green's functions of the dots electrons, here we show that 
those Green's functions can be computed self-consistently 
through a ^-derivable approximation, thereby ensuring 
that all conservation laws are respected. The method re- 
lies on a perturbation theory as formulated by Kadanoff 
and Baym^^, extended to the usual Keldysh contour to 
allow for a non-equilibrium situation, and will be referred 
to as the Kadanoff-Baym-Keldysh (KBK) approach. 

When applying this KBK formalism to the double-dot 
Cooper-pair beam splitter setup, we show that even weak 
interactions have an important effect on both current and 
cross-correlations in the CAR-dominated regime, and 
that there exists a competition between Coulomb inter- 
action and direct inter-dot tunneling. 

The outline of the paper is as follows. In section [TTl we 
discuss the general model for a set of quantum dots con- 
nected to multiple leads, and derive the corresponding 
expression for currents and current-current correlations. 
We then show in section IIIII how the KBK approach al- 
lows us to compute the single- and two-particle Green's 
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function. In section ITVl we apply this formalism to the 
case of the Cooper-pair beam-splitter. We conclude in 
section Ivl 



II. MODEL 



This averaging over the leads degrees of freedom allows 
us to express the relevant physical quantities in terms of 
the tunneling self-energy as well as the dots electrons 
single- and two-particle Green's functions. These are de- 
fined in Nambu-dot-Keldysh (NDK) space respectively 
as 



Our starting point is a set of quantum dots (labeled 
a), with energies e^, tunnel-coupled to multiple leads (la- 
beled j) with tunneling amplitudes taj. The leads have 
identical chemical potentials — fi and are character- 
ized by their voltage bias Vj and superconducting order 
parameter Aj (normal leads correspond to Aj = 0). The 
Hamiltonian for such a system is of the form 



(7) 



and 



(^1,^2,^3,^4) 



(1) 
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In terms of Nambu spinors dj^ — {(f^^ dai), the dots 
Hamiltonians read 



- (TKd.,.,(tr)rf„...(tr)'^i3.3(ir)4...(tr)), («) 

where Tk corresponds to the time-ordering operator 
along the Keldysh contour. 

In particular, the average current from the dot a into 
the lead j is readily obtained using standard methods^ 
as a function of the single-particle Green's function G 



{^^diaj^ + h.c.) + Uo^n^tnau =f TrNDK 



Ha = <^adtcr^Ja + ^ ( -^d'^a'^zAp + h.C 

(2) 

where we also introduced the inter-dot tunneling ampli- 
tude tap and the local Coulomb repulsion Ua- 
Similarly, the leads Hamiltonians are given by 



(r^ (X)Xq (g) (Tz) 



1- {^k(Tz + A,a^)m,k 



(3) 



where A. = ^7— 



-/i and we used the Nambu spinors ^'L 



(^'^j,^ "^j-ki)- Here and throughout the remainder of the 
text, we work with physical dimensions corresponding to 
h=l. 

We chose to transfer the voltage dependence onto the 
tunneling term using the Peierls substitution: Tja(t) 



i(7z Vn t 



The tunneling Hamiltonian is thus given by 



Hja^Y.{^krjait)da+h.C.). (4) 

fc 



Since the total Hamiltonian is quadratic in the leads 
electrons operators, we can integrate out these degrees 
of freedom. The effect of the leads is then captured by a 
tunneling self-energy t.T{ti,t2) = Y.j ^j(ii,t2). 



^J:A'lM2(^l,*2) - Tjaiitl) g]{tl,t2) 7jQ2(i2) 



mm 



(5) 



which depends on the tunneling amplitude written in 
Nambu-Keldysh space Tja{t) = ® Tja{t) as well as 
the standard non-interacting Green's function for lead 
electrons, given in Nambu-Keldysh space by 

g]:::-.,itut,) = ~^Y.{TK^,k.AtT)^k.JtT))■ (6) 

k 



For convenience, we also introduced the more compact 
notation /z^ = {ai,ai,r]i}, condensing the Nambu-dot- 
Keldysh components into a single index. 



X / dt' {G{t,t')tj{t',t)~tj{t,t')Git',t)) 



(9) 

where the matrix la is defined in dot space, with ele- 
ments [Ia]aia2 ^ ^aaiSaas, and TrNDK corresponds to 
the trace in the full Nambu-dot-Keldysh space. 

Similarly the current correlations can be expressed in 
terms of the two-particle Green's function K as 



dtidt2 ^ 



MlA'2 



~ Si.f/Ji (i, ^l)^j,/i2i'' (^2 , ^ )K ^-^^1,1 ^^^{tijt ,^,^2) 
~ t)'^j,i^'p,2 (^', ^2)-^i//^2A'li'' (^, ^2, ^1, 

{t,t',ti,t2)\, (10) 

with the labels v = {a, a,!]} and i^' = {cr', /3, 77'}. 



III. KADANOFF-BAYM-KELDYSH 
FORMALISM 

In order to proceed further, we need to evaluate the 
fully interacting dots electrons Green's functions G and 
K. This is achieved through a perturbative expansion 
in the Coulomb interaction Ua, within a so-called $- 
derivable approximation^^. Such a conserving approxi- 
mation relies on the truncation of the Luttinger- Ward- 
functional $ [G] . The latter is a function of the fully 
interacting single-particle Green's function G, and cor- 
responds to the sum of all closed-loop two-particle ir- 
reducible diagrams. Unfortunately, there is no simple 
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FIG. 1: Diagrammatic representation of (a) the contri- 
butions to the Luttinger-Ward functional, (b) the Bethe- 
Salpeter equation in the horizontal particle-hole channel, and 
(c) the two-particle Green's function in terms of the single- 
particle one and the vertex function. The full lines correspond 
to full single-particle Green's functions for the dot electrons, 
the disks represent bare vertices. The large circles and squares 
correspond to the irreducible and full vertices respectively, 
while the shaded block stands for the two-particle Green's 
function K defined in the text. 



closed form of the Luttinger-Ward functional, which is 
usually represented in a more convenient diagrammatic 
form (see Fig. [1^). 

The first striking feature of the Luttinger-Ward func- 
tional is that it serves as a generating functional for 
the interacting self-energy. Indeed, taking the func- 
tional derivative of $ [G] with respect to the one-particle 



Green's function G leads to 



(5$ [G] 



(11) 



Note that the resulting expression for the self-energy is 
itself a function of the full propagator, so that G has to 
be derived self-consistently. This is achieved by solving 
the Dyson equation 

G{h,t2)=G°ih,t2) + J dt^dU&ihM) 

X [tc{h,U) + trihM)] G{U, t2), (12) 

where & is the dots electrons Green's function in absence 
of both interactions and tunneling to the leads. 

Further differentiating the Luttinger-Ward functional 
gives access to more involved many-body objects. In par- 
ticular, one can derive an expression for the irreducible 
vertex from a second order functional differentiation. In 
the horizontal particle-hole^i channel, it reads 



A„ 



i(il, ^2, ^3, ii) 



J2$ [G] 



(13) 



This in turn allows us to write the full vertex function F 
by self-consistently solving the Bethe-Salpeter equation 
in the corresponding irreducibility channel (see Fig. [TJ)) 



r/iiAi2Ai3A'4(^lj^2,i3,i4) — A^j^2P3M4 (^17^2,^37^4) 

+ j dt^.-.dts ^ A^jp2pg^g(ti,<2,i5,^6)G^5M7(^5,i7)G^8M6(*8,^6)f^7A'sM3A'4(*7,i8,i3,i4)- 
M5---MS 

(14) 

The two-particle dots electrons Green's function K is then readily obtained from 

-^A'lA'2P3A'4 (^1,^2,^3,^4) = ^^1^4(^1,^4)0^2^3(^2,^3) — G ^-^^^{ti,t^)G ^2^^{t2,t4) 

+ ^ j (^^5 ■ • ■^^i8G^ip5(ii,i5)Gp2M8(^2,i8)fp5^g^7pg(i5,^6,^7,^8)Gp7M3(^7,i3)Gpg^t4(t6,^4)- 

(15) 

I 

It follows that from a given truncation of $ [G] , one can only one-particle self-consistent since the Luttinger-Ward 
determine both the full single- and two-particle Green's functional is a function of G only, 
function, respectively via the derivative of the Luttinger- 
Ward functional and the self-consistent solution of the 
Dyson equation, then through the second derivative of 
<& [G] , after solving the Bethe-Salpeter equation. Note 
that although both Dyson and Bethe-Salpeter equations 
have to be solved self-consistently, the proposed scheme is 
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FIG. 2: Cooper-pair beam-splitter setup: (a) representation 
of the experimental device where a nanowire/nanotube is cou- 
pled to normal and superconducting leads, and (b) schematic 
representation of the double-dot model considered in the text 
(with the corresponding notations). 



IV. APPLICATION TO THE COOPER-PAIR 
BEAM-SPLITTER 

Let us now illustrate this procedure on the particu- 
lar example of the double-dot Cooper-pair beam splitter, 
whose setup is recalled in Fig. O 

A. Derivation of the self-energy and vertex 
function 



so that our conserving approximation amounts to writing 

^ J . <j<j' aa' 

(16) 

where a — —a. 

Following the procedure outlined in the previous sec- 
tion, we first obtain the interacting self-energy, which in 
frequency space writes 

(17) 

One readily sees from this expression that at this level 
of approximation the interacting self-energy in Fourier 
space is frequency independent, and can thus be viewed 
as a renormalization of the dots energies. However, 
this constant still needs to be evaluated self-consistently, 
through the Dyson equation (|12p . as it depends on the 
full single-particle Green's function G. 

The irreducible vertex is further obtained following 
Eq. (1131) . and takes the following simple frequency in- 
dependent form in Fourier space 



Our system now consists of a superconducting lead at 
potential V5 = 0, coupled to two quantum dots (labeled 
1 and 2), which are each connected to a different normal 
electrode (at potential Vl and Vn respectively). We will 
be interested in the branching currents Ili and 7^2 as 
well as their cross-correlations. Here Ili corresponds to 
the current between dot 1 and the left normal lead, while 
Iji2 flows between dot 2 and the right normal lead. 

In what follows we focus on the first term in the di- 
agrammatic expansion of the Luttinger-Ward functional 



A. 



^j^mmx S U S S S 



X a 



(18) 



where we used the energy conservation to introduce 
A(cji,W2,W3,aj4) — 2ttS {u!i — UJ2 — UI3 + UJ4) A. The re- 
sult of Eq. ([T^ is nothing but the bare interaction ver- 
tex where one recognizes the properties of locality in time 
(all ti and rji are equal) and space (all ai are equal), as 
well as the spin conservation (the sum of incoming and 
outgoing spin projections is zero). 



J 



Substituting this expression for the irreducible vertex back into the Bethe-Salpeter equation (jl4[) and moving to 
Fourier space, one notices that the frequency-dependent full vertex function depends on a single frequency argument, 
so that by further introducing r(a;i, a;2, W3, 1^4) — iirS {coi — a;2 — + W4) F {uj2 — t^i), we are left with 



^ E 



(CJ2-W1). (19) 



In order to extract an explicit expression for the full 
vertex, one needs to invert Eq. (dH), which cannot be 
straightforwardly achieved since we are dealing with 4*^ 
order tensors. To circumvent this issue, we propose 
to construct a matrix representation of the full and 
irreducible vertices in an enlarged NDK(g)NDK space, 
by combining two by two the four labels of these ten- 
sor elementsi^ Indeed, the matrix representation of the 
single-particle Green's function relies on identifying a sin- 
gle coordinate in NDK space ni (running from 1 through 



8) out of the values of the Nambu, dot and Keldysh com- 
ponents {(Ti, Qfi, ?7i}. In the same spirit, we construct en- 
larged 8^ X 8^ matrices to represent the full and irre- 
ducible vertices, in which the coordinates are obtained 
from a linear combination of the rii according to 

^-fQiQ2 — L Jmn, -Laia2a3a4 — nin2n3n4 — I i Nt No 

(71(72 ^ ^ (Tl (72 (73(74 ^ ^ 

(20) 

with Ni = 8(ni - 1) -|- n2 and iV2 = 8(n3 - 1) -I- n^. The 
resulting 64 x 64 matrices, which we hereby denote with 
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FIG. 3: Currents Il\ and I fa (in units of Jo = eA/ft) as a function of Vl, for different values of the interaction I] = 
0, 0.05, 0.1, 0.15A, and the parameters (in units of A) /3 = 100, ei = 0.5, £2 = -0.5, Vr = -0.7 and t^x = tsx = ts2 = = 0.2. 
The upper panels correspond to td ~ 0, while the lower panels are computed for td ~ 0.2A. 



brackets, ultimately lead to the following explicit form 
for the full-vertex function 

[r{uj2-cji)] = [A] {[i]-[n(c^2-c^i)] X [A\y\ 

(21) 

where we introduced the generalized polarization bubble 
[n] as 

moj)] =J^ [cm G^{n + c.)] , (22) 
and is the transpose of the matrix G in NDK space. 



B. Numerical results 

We computed numerically the currents and crossed 
correlations as a function of the left lead voltage Vl (with 
fixed Vr), for different values of the interaction param- 
eter Ui = U2 = U. We focused on the most favorable 
regime for CAR, i.e. when the dots energies are chosen 
antisymmetric with respect to the superconductor chemi- 
cal potential, and considered the situations of a negligible 
and a strong direct tunneling. 

We first solve self-consistently the Dyson equation in 
order to compute the single-particle Green's function G 
and derive the currents Ili and 1^2 from Eq. ([9]), which 
we present in Fig. [3] In absence of direct tunneling, the 



two currents are almost identical, while for a finite t^, 
the current difference is significant, due to the opening 
of a new conduction channel between the normal leads. 

As a function of U, one readily sees that the amplitude 
of all currents gets reduced, already for small values of the 
interaction strength. This reduction does not equally af- 
fect the whole range of voltages. In particular the region 
of high voltages, where the current is dominated by elec- 
tron co-tunneling is only marginally modified. However, 
the current reduction in the presence of interactions is 
sizable in the low-voltage CAR-dominated regime. More- 
over, the effects are qualitatively similar whether a direct 
tunneling is present or not. These are purely static single- 
particle effects, and can be understood as a detuning of 
the resonances in the dots density of states, ^"^ a direct 
consequence of the constant interacting self-energy. Note 
that while the latter is frequency- independent, it does 
depend on Vl (through self-consistency), explaining why 
different regions of the curves are not equally affected. 

The most interesting and non-trivial results come from 
the calculation of the current cross-correlations as both 
one- and two-particle effects are present in this case. In 
what follows, we focus on the zero- frequency current- 
current correlation 5^^73(0) defined as 



5,„,,73(0) - / dt ( (/^(i)/+ (0)) - {hUt)){iA0)) 



(23) 



I 

Substituting the expressions ^ for the branching currents and (|10p for the current-current correlators, and using 
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Eq. (jl5|) to replace the two-particle Green's function in terms of the full vertex, the zero-frequency correlations become 



+ 



/diiJi 



duj2 



^15, Me, Mr, Ms 



(Si^iy^jG^l^jG^^y) (Wl) (Sj;^2l^'G'iy'^sG'pg^2) (^2) — (Si^p^^Gi/p^G^^^i) (Wl) (Sj,iy>2G'M2A'8G'/i6l'') ('^2) 

+ (Si,^iiyG^^5G^7^i) (wi) (Sj,^2^/Gi,/pgG^|,^2) (a;2)] |. (24) 

I 



Using the matrix representation introduced above, the 
cross-correlations can be recast in a much more manage- 
able form as 



2^TrNDK [^LiG^R2G] i^) 
'/ ^^TrNDK^{[f(c.2-c.i)] 
(GS^iG)(^i)®(GS+2G)'^M]}, (25) 

where we introduced the trace Tr^Qj^^ over our enlarged 
NDK(g)NDK matrices and defined the compact notations 



TTai = i (1 ± T^) ®r 



(26) 
(27) 



The first term on the r.h.s of Eq. is similar to 

the expression obtained in the absence of interactions j^'^ 
only here it involves the full dots electrons single-particle 
Green's function. The second term however corresponds 
to truly two-particle effects. 

In Fig. m we show the results obtained for the zero- 
frequency current-current correlations in the absence 



of direct tunneling. Like in the non-interacting case, 
one can typically isolate two relevant regimes of volt- 
ages. At low voltage (Vl < ei), transport is domi- 
nated by CAR processes and we observe positive cross- 
correlations, which are approximately constant. Increas- 
ing the interaction parameter leads to strongly enhanced 
correlations in this region. This enhancement cannot be 
solely attributed to the renormalization of G due to in- 
teractions. We could check that the vertex-dependent 
term in Eq. ((25|) contributes substantially to the final re- 
sult, therefore suggesting that interaction- induced two- 
particle effects are particularly important in the Cooper- 
pair beam-splitter setup. 

In the region of high voltage {Vl > ei), we came across 
a somewhat counter-intuitive result. In this regime, 
CAR is greatly unfavored and one expects electron co- 
tunneling to be the dominant process, as confirmed by 
the current plots of Fig. [3] where Ili and 1^2 are opposite 
in this region therefore signaling the transfer of electrons 
from the left to the right normal lead. One would thus 
expect to observe negative correlations, as was the case 
in the absence of interactions. However, paying close at- 





FIG. 4: Zero-frequency current cross-correlations (in units of 
So = e'^A/h) as a function of Vl, for U = 0, 0.05, 0.1, 0.15A, 
and the same parameters as in Fig. [3l in the absence of direct 
tunneling. 



FIG. 5: Zero-frequency current cross-correlations (in units of 
So = e^A/h) as a function of Vl, for U = 0, 0.05, 0.1, 0.15A, 
and the same parameters as in Fig. |3l in the presence of direct 
tunneling, td = 0.2A. 
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tention to the high- voltage region (see the magnified area 
of Fig. S]), we see that the effect of interactions, although 
weak, is sufficient to flip the overall sign, resulting in 
slightly positive correlations. 

When a direct tunneling is allowed between dots, there 
is still an important enhancement of positive correlations 
at low voltage, even for weak interaction strength (see 
Fig.[5])- Also, in contrast with the constant behavior ob- 
served for td ~ 0, SL more pronounced feature develops 
around — — Ci. This is reminiscent of what was ob- 
served in the non-interacting case, and can be attributed 
to density of states effects^^. For voltages Vl > ei, inter- 
actions have a weaker effect, and tend to reduce the am- 
plitude of cross-correlations, without causing any change 
of sign. 

Comparing the results of Figs.|3]and[Sl it becomes clear 
that while a finite td tends to shift the correlations to- 
wards negative values, a finite U have the opposite effect 
of shifting them towards positive ones, thus revealing the 
competition between interactions and direct tunneling. 

V. CONCLUSION AND OUTLOOK 

In conclusion, we presented a conserving approach 
to study currents and their correlations in interacting 
out-of-equilibrium mesoscopic systems. This Kadanoff- 
Baym-Keldysh formalism amounts to computing the 
single- and two-particle dots electrons Green's functions 



from a given truncation of the Luttinger-Ward func- 
tional. These are then used to express the branch- 
ing currents (in the same spirit as the Meir-Wingreen 
formula^^) as well as their correlations, therefore extend- 
ing the Fisher-Lee/Landauer-Biittiker formula^- for the 
noise to the case of an interacting system. Applying this 
formalism to the double-dot Cooper-pair beam splitter 
setup, we could study the effect of Coulomb interaction 
on the currents and cross-correlations, showing that even 
weak interactions have an important effect in the CAR- 
dominated regime, in particular favoring positive cross- 
correlations in the antisymmetric case. Our results also 
exhibit the competition between local Coulomb interac- 
tion and direct inter-dot tunneling, which affect the cur- 
rent correlations in an opposite manner. 

The present work can be extended in various ways. 
For example, the formalism can also treat a non-local 
Coulomb interaction between dots. This would yield new 
terms in $ [G] , and therefore in the self-energy and ir- 
reducible vertex, which can be easily accounted for. But 
the most natural extension consists in including a second 
set of diagrams into the Luttinger-Ward functional (see 
Fig. [1^). In particular, this second order contribution 
would lead to a newly acquired frequency dependence 
of the interacting self-energy and irreducible vertex, and 
thus to new features associated with these retardation ef- 
fects. However, this would also make the Bethe-Salpeter 
equation no longer invertible analytically, therefore re- 
quiring a self-consistent numerical solution. 
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